Absorbing states and elastic interfaces in random media: 
two equivalent descriptions of self-organized criticality 
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We elucidate a long-standing puzzle about the non-equilibrium universality classes describing 
self-organized criticality in sandpile models. We show that depinning transitions of linear interfaces 
in random media and absorbing phase transitions (with a conserved non-diffusive field) are two 
equivalent languages to describe sandpile criticality. This is so despite the fact that local roughening 
properties can be radically different in the two pictures, as explained here. Experimental implications 
of our work as well as promising paths for future theoretical investigations are also discussed. 

PACS numbers: 05.50.+q,02.50.-r,64.60.Ht,05.70.Ln 

The concept of self-organized criticality (SOC) has 
been proposed to account for the emergence of scale in- 
variance in Nature Its main tenet is that in the 

presence of slow driving and fast dissipation, acting 
at infinitely-separated time scales, many systems self- 
organize (without any explicit tuning of parameters) to 
a critical state Some archetypical examples of SOC 
are provided by sandpile toy models, in which grains are 
slowly added, locally redistributed on a fast timescale 
whenever an instability threshold is overcome (generating 
avalanches of toppling events), and finally dissipated at 
the open boundaries. Upon iteration, this process leads 
to a critical steady state. 

SOC can be related to standard (non-equilibrium) crit- 
icalphenomena by defining the "fixed energy ensemble" 
[H, 3 m which driving and dissipation are switched off, so 
that the number of grains (or "energy") is conserved. Us- 
ing this quantity as a control parameter, a standard (i.e. 
non self-organized) phase transition is observed: for large 
energy densities, there is a finite density of active (top- 
pling) sites, whereas the system ends in a frozen (stable) 
state at low densities. It has been shown that the critical 
point separating these two regimes occurs at the value 
of the energy density at which the system self-organizes 
when subjected to slow-driving and boundary dissipa- 
tion 0, Subsequent debate has attempted to eluci- 
date which non-equilibrium universality class stochastic 
sandpiles belong to. As detailed below, two alternative 
solutions have been proposed. 

Sandpiles were first related to interfaces in random me- 
dia |6j. In this language, the interface height, h(x,t) is 
the number of times a given site x has toppled up to time 
t, and frozen states correspond to pinned interfaces. The 
resulting pinning-depinning transition was argued to fall 
in the quenched Edwards- Wilkinson or linear interface 
model (LIM) class, described by @,[1]: 



This correspondence was recently proven exact between 
one particular sandpile model and one member of the 
LIM class [l(J but, in general, it is only approximate, as 
some noise-correlations need to be neglected to establish 
a full correspondence with Eq.JT]). 

Alternatively, sandpile models have been rationalized 
as systems exhibiting an absorbing-state phase transition 
Indeed, in the fixed energy ensemble, a stable config- 
uration is one of the infinitely-many absorbing states in 
which the system can be trapped forever, whereas activ- 
ity never ceases above the critical point. The correspond- 
ing universality class is not the prominent directed per- 
colation (DP) class, but is characterized by the coupling 
of a DP-like activity field to a conserved, non-diffusive, 
auxiliary field (the "energy") dEl. often called C - DP 
(or also Manna [13 ]) class, it is characterized by the fol- 
lowing set of Langevin equations: 
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d t h(x, t) = V 2 h(x, t)+F + r)(x, h) 



(1) 



where F is a force and rj{x, h) is a quenched white noise. 



where p is the activity field, <j) the background energy 
field, and r)(x,t) a Gaussian white noise [4, lllj. 

The validity of both of these alternative pictures has 
been (partially) backed by numerical measurements of 
critical exponents but, in general, they have not been 
proven to be correct so far. But, if both of the pic- 
tures are right, a remarkable consequence follows: depin- 
ning transition of LIM-class interfaces in random media 
and the C-DP class absorbing phase transition should be 
equivalent, even if they look rather different (for instance, 
one involves quenched disorder and the other does not). 
Here, we explore this issue and the more general question 
of whether any depinning interface universality class has 
an equivalent absorbing phase transition class. 

Only a few works have approached the connections be- 
tween these two pictures. In 13|, interfaces were con- 
structed from a DP class model, using the cumulated lo- 
cal activity as the interface height. Anomalously-rough 
interfaces, characterized by a positive local-slope expo- 
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TABLE I: Some of the measured critical exponents of the C- 
DP/LIM class vs space dimension d @, B EE E3, ■ The A- 
and B-scaling values for k are from our own present simula- 
tions. 
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nent k defined by ((V/i) 2 ) ~ t 2fi [lj], were found at criti- 
cality. These anomalous interfaces are not related to any 
known interface class. Focusing on SOC sandpiles, Alava 
and Munoz [l5| argued heuristically that the LIM and 
C-DP classes could be identified with each other (using 
also h(x,t) — Lp(x,s)ds to relate Eq.flT} to Eq.@) al- 
though a one-to-one mapping could not be rigorously es- 
tablished. However, this conclusion was later challenged 
by Kockelkoren and Chate (KC) [l6| who found that the 
k exponent takes completely different values for LIM in- 
terfaces and for interfaces constructed from models in the 
C-DP class. In particular, the constructed interfaces are 
anomalously-rough below the upper critical dimension, 
i.e. for space dimensions d < 4, while LIM interfaces have 
K < (not anomalous) for d > 2. In d = 1, both types 
of interfaces are anomalously rough, but in a manifestly 
distinct manner, i.e. different values of k (Fig. [1]). This 
led KC to conclude that LIM and C-DP classes cannot 
be equivalent, even if all the other recorded "standard" 
critical exponents (as 9 and z, see definitions below and 
table I, and others) take "almost indistinguishable val- 
ues" in these two problems, (which could, in principle, 
be attributed to a numerical coincidence [161].) 

While the discrepancy in values of k is unquestionable, 
simulating directly Eqs.© (using the method in [l7| ) we 
find all the other C-DP exponents to be indistinguishable 
from their counterparts in the LIM class (Table [TJ as well 
as from the corresponding values in stochastic sandpiles. 

In this Letter, we show that depinning transitions of 
LIM interfaces and C-DP absorbing phase transitions are 
indeed two equivalent descriptions of SOC sandpiles in 
spite of the discrepancies in k- values, which we explain. 
We show, using a combination of numerical results and 
scaling arguments, that there is a unique universality 
class and that differences in k- values stem from diverging 
local fluctuations, inherent to the absorbing state picture, 
which do not affect other long-distance properties. 

Let us start by clarifying the origin of the two pos- 
sible values of k and the scaling laws they obey by us- 
ing simple scaling arguments. First, since h — J dtp(t) 
its scaling dimension is [h] ~ t 1_e , where 9 is the den- 
sity (or interface- velocity, recalling that dth = p) crit- 
ical time-decay exponent: (p(t)) ~ t~ e . This leads to 
[Vft] ~ t -i/z+x-e^ where [ V -i] _ t i/z defines the dy _ 
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FIG. 1: Typical one-dimensional interfaces for a system 
of 2 15 sites at t = 10 5 from flat initial conditions. Top: 
Leschhorn automaton (LIM class, k « 0.17(1)). Bottom: in- 
terface constructed from the Manna sandpile |12|] (C-DP class, 
k » 0.43(1)). 



namical critical exponent z, and therefore 
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("A scaling" from now on). (3) 



For d = 1, plugging the values 9 « 0.13 and z s» 1.42 of 
the LIM or C-DP class 0, H, [3, [3 into this expression 
leads to k sa 0.17 which is indeed the value measured for 
LIM class interfaces 8]. In higher dimensions, this scal- 
ing law yields zero (with possible logarithmic corrections 
in d = 2) or negative k values, i.e. no anomalous scaling, 
as indeed observed in simulations [lq |. 

On the other hand, assuming that interface heights at 
adjacent sites are asymptotically uncorrelated (this will 
be justified after) we have V/i ~ Vh ~ t( 1_e )/ 2 and hence 
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k b = 



("B scaling") 



(4) 



KC observed that B-scaling is verified by many interfaces 
constructed from microscopic models at absorbing phase 
transitions (l6j . Our own simulations (not shown) ex- 
tend this result to different sandpile models (simulated 
in the fixed e nerg y ensemble): the constructed interfaces 
of the Manna [l2j , Oslo @| , and Mohanty-Dhar mod- 
els show B-scaling at criticality. 

To shed some light on the physical reason for the 
existence of two different k values, let us consider the 
Leschhorn automaton, a LIM-class model showing A- 
scaling. It is a discretization of Eq.{T]): an integer- valued 
height advances at each site x following: 

h(x) -> h(x) + 1 iff V 2 h + F + <q{x, h{x)) > 0, (5) 

where the (discrctized) Laplacian is computed using the 
nearest- neighbors of a;, and rj = ±1 with respective prob- 
abilities p and 1 — p. Consider now the Manna sandpile, 
a C-DP class model whose local rule is: if two or more 
grains are present at a given site, distribute two of them 
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FIG. 2: Time series of interface squared-gradient and activity- 
density (or velocity) at criticality. (a) Eqs.([2} using h, (A- 
scaling, k ~ 0.17(1)) or h (B-scaling, k m 0.43(1) ). Parame- 
ters: L = 2 15 , b = w = a 2 /2 = 1, D = D E = 0.25, time-mesh 
0.1; critical point a c = 0.86452(5). (b) Leschhorn automaton 
with L = 2 15 , F = 0, at the critical point p = 0.80078(5), and 
using h as in Eq.JT} (A-scaling, k » 0.17(1)) or h as in Eq.© 
(B-scaling, n » 0.43(1)) with p c = 0.76935(5). 

randomly to the nearest-neighbors [1, [Hj]. In this case, 
the interface h(x) encoding the number of times a site 
has toppled since t = shows B-scaling, as said before. 
It can be can be expressed in terms of Zi n (x) and z ut(x), 
the cumulated number of particles respectively received 
from and given to the nearest-neighbors. For example, 
in d — 1, z out (x) — 2h(x) while z- ln can be expressed as 
the sum of a mean flux h(x + 1) + h(x — 1) plus a fluc- 
tuating part, £(x,h(x)) : indicating stochastic deviations 
from this mean [6|, Il5l |l6| . The toppling condition can 
be written as zq(x) + Zin(x) — z ou t(x) > 2 (where zq (x) is 
the initial number of grains) and thus , expressed in terms 
of the following advancement rule: 

h(x)^h(x) + l iff V 2 h- l + z (x) +£(x,h) > . (6) 

This is very similar to Eq.© but, as noticed in [l6|, there 
is a crucial difference: whereas in Eq. (|5|) rj(x,h(x)) is a 
bounded, dichotomous, delta-correlated noise, the noise 
term £(x,h(x)) in Eq. ^ is a sum of random variables 
(a unit is added or subtracted for each toppling) whose 
amplitude, by virtue of the central limit theorem, behaves 
like the square root of the average of h(x + 1) + h(x — 1), 
and is therefore diverging in time: (£,(h) 2 ) ~ t 1 ^ 8 . In 
turn, this divergence has to be compensated by the fluc- 
tuations of the Laplacian term in Eq. ([6]) (since the term 
Zq {x) representing the initial condition should be irrel- 
evant in the long-time limit and is anyhow bounded). 
This is at the origin of the strong fluctuations present 
in the constructed interface (B scaling) but absent in 
the LIM class (see Fig[T]). At this point, one clearly 
appreciates the qualitative difference between LIM and 
C-DP-constructed interfaces, which occurs despite of the 
fact that both classes share numerically-indistinguishable 
(standard) critical exponents, as said before. We now 
show that this is not the end of the story, and that we can 



construct A-scaling interfaces from C-DP class models, as 
well as modify LIM-class models to obtain B-scaling. 

Our first evidence showing that A-scaling and B- 
scaling can both be compatible with a unique uni- 
versality class was provided by numerical integrations 
of Eqs.@. Constructing an interface, as before, via 
h(x,t) — J dsp(x,s) where p is the continuous activ- 
ity field, we obtain clear A-scaling (Fig. [2j, in contrast 
with the B-scaling heretofore always observed with mi- 
croscopic models. Next, mimicking microscopic models 
in which the interface advances by one unit whenever 
a site is active, we constructed a different interface for 
Eqs.© through h(x,t) — J Q dsQ(p(x, s)) where O is the 
Heaviside step function; this new interface advances by 
one unit whenever there is some non-zero activity, re- 
gardless of its magnitude. Strikingly, B-scaling is then 
observed (Fig. [5]). Thus both A- and B-scaling interfaces 
have been constructed at the same absorbing phase tran- 
sition point: they correspond to slightly different observ- 
ables. We have reached a similar conclusion for the Oslo 
sandpile model by taking advantage of a recent result 
by Pruessner [lfj who constructed an exact mapping be- 
tween this particular sandpile and the Leschhorn automa- 
ton. The local rule for d — 1 is as follows: distribute one 
grain to each nearest-neighbor whenever the local height 
threshold is passed, this local threshold being randomly 
reset to be 1 or 2 grains after each toppling. To achieve 
an exact mapping, Pruessner showed that it is crucial to 
use the more symmetrical fv(x) — h(x + l) + h(x — l), i.e. 
the accumulated number of times a given site has been 
charged by its neighbors rather than the accumulated 
activity at the site itself (the definition of h). Indeed, 
using tv(x) for the Oslo model eliminates the diverging 
noise in Eq.((6]) and yields A-scaling in numerical simu- 
lations, whereas using h(x) we observe B-scaling (results 
not shown). We have been able to extend easily this 
procedure to other sandpiles as the Mohanty-Dhar one 
[19(. For other sandpiles with less symmetric redistribu- 
tion rules as, for instance, the Manna one [13] this can 
be much more complicated. In this last, the 2 toppling 
grains at any site can go to the same neighbor. This in- 
troduces an extra noise that needs to be subtracted by an 
appropriate (and intricate) definition of the height vari- 
able (different from h and fv) to get rid of intrinsic local 
fluctuations and disentangle the hidden A-scaling. These 
results show that a unique universality class is compati- 
ble with different values of k, depending on microscopic 
details and/or the definition of the height variable; dif- 
ferent K exponents correspond to different, though very 
similar, observables. Note also that the same definition 
of the interface can lead to the two different types of 
scaling depending on the rules of sandpile under study. 

To close the loop, we now show that for standard inter- 
face depinning transitions it is also possible to generate 
two different k values without affecting other exponents. 
Let us define dth(x,t) = dth(x,t)(l + a{x,h)), where 



4 



the quenched noise a is or 1 with probabilities p and 
1 — p, respectively, so every time the original interface 
advances, a noise variable is added to h. Both interfaces 
are related by h(x, t) = h(x, t) + a(x, h) where a(x, h) is 
an accumulated noise summing up all values of a(x, h) 
at x up to height h. In terms of h(x,t), Eq.([T]) becomes 

d t h = [S7 2 (h(x,t)-5{x,h)) + r]{x,h)} x [l + a(x,h)]. (7) 

Simulating Eq.0 we observe B-scaling for the h- 
interface, while removing the cr-noise we readily recover 
A-scaling for the ft,-interface (Fig. [2|) , while all the other 
exponents coincide for both interfaces. Naive power- 
counting for Eq.([7]) shows that V 2 a(x, h) is an irrelevant 
higher-order noise, and so is the term a(x, h) added to 
1 [tJ- Hence, upon coarse-graining, Eq. (|7j) flows towards 
the standard LIM renormalization group fixed point Q, 
justifying that all universal critical exponents should co- 
incide in Eq.|T]) and Eq.([7]) in accordance with our nu- 
merical results. The inclusion of the extra noise (a 
higher-order correction to scaling) is thus able to alter 
the value of k, intensifying anomalous behavior, but not 
standard long-distance critical exponents, which are con- 
trolled by a unique renormalization group fixed point. 
The two different values of k correspond to two different 
height variables, h and h, differing by a diverging noise. 
A full understanding, within the renormalization group 
perspective, of how local anomalous roughening proper- 
ties are affected by an otherwise irrelevant noise remains 
a challenging task. 

The fact that a given universality class can be compat- 
ible with different types of local roughening (amenable to 
experimental analysis) is also of interest in general stud- 
ies of fluctuating interfaces [14j. Indeed, one can show 
that x an d Xioc being respectively the global and lo- 
cal saturation roughness exponents are related by \ — 
Xioc + zk 14], which we have numerically verified for 
all models studied here. Given that our results indi- 
cate that local roughening properties, as encoded by k 
( or Xioc), can adopt different values, depending on the 
presence or absence of local fluctuations, while X: as an 
other non-local exponents, is universal, it is intriguing 
that experiments on the propagation of fracture cracks 
in wood seem to lead to the opposite conclusion. 
See also the nice recent experiments measuring the ex- 
ponents reported here for self-organized superconductors 
[2l| . More experimental work along these lines would 
be most welcome, since experimental realizations of ab- 
sorbing phase transitions (even for the paradigmatic DP 
class) are barely existing [221 ]. 

On the theoretical side, since the C-DP class and the 
LIM class are two faces of the same problem, both share 
the same upper critical dimension d u = 4. This re- 
sult is in contradiction with the perturbative approach 
for Eq.([2|) in [23[. In fact, field-theoretical treatments 
of the LIM class demand a functional renormalization 



group calculation and one finds that the correlator 
of the quenched noise develops a cusp in the h- variable. 
Given our results, it would be very interesting to sort 
out the analogue of all this in the absorbing phase tran- 
sition picture, as well as attempting a non-perturbative 
renormalization group approach for such a case. 

In summary, depinning transitions of linear interfaces 
in random media and absorbing phase transitions in the 
C-DP class are two equivalent descriptions of sandpile 
self-organized critical points. This clarifies the issue of 
universality in stochastic sandpiles and the connection of 
SOC to standard non-equilibrium phase transitions and 
opens the door to new and exciting research lines. 
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